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Free-Free  Transitions  in  Collisions  between  Slow  Electrons 
and  Neutral  Oxygen  Atoms 


Reprinted  (with  minor  modifications)  with  permission  from 
MFree-free  transitions  between  slow  electrons  and  neutral 
oxygen  atoms11  by  Sunggi  Chung  and  Chun  C.  Lin,  Physical 
Review  A,  Volume  51,  p.1221  -  p.1229,  February  1995. 
Copyright  1995  The  American  Physical  Society. 
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1.  INTRODUCTION 


During  an  encounter  with  an  atom,  an  electron  may  emit  (Bremsstrahlung)  or  absorb  (in¬ 
verse  Bremsstrahlung)  a  photon  in  a  manner  that  the  atom  is  not  affected  energetically.  Such 
free-free  transitions  of  electrons  are  important  source  of  emission  and  absorption  processes 
in  the  stellar  atomospheres  and  the  low-temperature  plasmas  [1-3]. 

The  free-free  transitions  with  the  neutral  atomic  oxygen  as  a  host  have  been  studied 
theoretically  by  several  groups  [4-7].  Since  the  free-free  transitions  assume  a  particular 
importance  at  low-incident  electron  energies,  usually  below  the  first  excitation  threshold, 
the  target  polarization  is  an  important  factor  to  consider  in  the  theoretical  calculation. 
Thus,  in  the  work  of  Mjolsness  and  Ruppel  [4]  and  of  Geltman  [5]  an  empirical  polarization 
potential  is  included  in  the  calculation  of  the  continuum  wave  functions.  John  [6]  obtained 
the  absorption  coefficients  from  the  momentum-change  cross  sections  by  means  of  a  formula 
suggested  by  Dalgarno  and  Lane  [8].  Later  John  and  Williams  [7]  adopted  the  multichannel 
theory  of  Morgan  and  John  [9]  to  calculate  the  absorption  coefficients,  utilizing  the  results 
of  a  three-state  close- coupling  calculation  of  the  electron-atom  scattering  [10]. 

In  this  work  we  account  for  the  target  polarization  by  means  of  the  method  of  polarized 
orbitals  which  Temkin  [11]  devised  to  study  the  e-0  scattering  and  later  extended  to  the 
e-H  problem  [12].  Henry  has  also  applied  this  method  to  study  the  elastic  scattering  of  the 
e-0  system  and  photo  detachment  from  0“  [13].  In  the  polarized-orbitals  method,  the  target 
polarization  is  made  a  part  of  the  atomic  wave  functions  which  depend  parametrically  on 
the  incident-electron  coordinates.  The  polarized  orbitals  in  turn  present  the  (polarization) 
potential  to  the  scattering  electron.  Thus,  we  are  in  a  position  to  compare  the  results  based 
on  three  different  ways  of  handling  the  target  polarization  -  by  the  empirical  polarization 
potential,  by  the  multi-state  close-coupling,  and  by  the  method  of  polarized  orbitals. 

In  Sec.  2,  we  briefly  discuss  the  basic  formulation  of  free-free  absorption  and  emission 
coefficients,  and  describe  the  method  by  which  the  free-electron  functions  are  obtained.  The 
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polarized-orbitals  method  of  Temkin,  from  which  the  polarization  potentials  are  derived,  is 
described  in  Sec.  3.  A  short  description  of  the  numerical  procedure  is  contained  in  Sec. 
4.  In  Sec.  5,  we  compare  our  results  with  the  previous  works  in  some  detail,  particularly 
with  the  work  of  John  and  Williams  [7],  since  our  approach  and  theirs  may  be  viewed  as 
complementary.  Section  6  concludes  this  report. 

2.  THEORY 

2.1  Absorption  and  Emission  Coefficients 

We  study  the  free-free  transition  process 

e-fy  +  0(2p4  3P)  -»  e-(kf)  +  0(2p4  3P )  ±  ftw,  (1) 


in  which  a  photon  is  emitted  ( +hu> )  or  absorbed  The  conservation  of  energy  gives  the 

relationship 


(kihy  _  {kfhy 


dh  tuv. 


(2) 


2m  2m 

It  is  appropriate  to  construct  the  wave  functions  of  the  (e“+0)  system  that  are  eigenfunctions 
of  the  toted  orbital  (L)  and  spin  (S)  angular  momenta  if  we  neglect  the  spin-orbit  interaction. 
Then  the  transitions  of  interest  are  of  the  type 


(2p4(SiZi)MSX]  -  [2p4(SiL1)kJ£SL'\t  (3) 


where  Si=Li=l  for  the  3P  state,  and  {1,1')  are  the  partial  waves  of  the  incident  and  scattered 
electrons.  The  dipole  matrix  elements  for  the  above  transition  have  been  worked  out  [16]. 
Using  these  results  the  formula  for  the  absorption  cross  sections  cah,[\,2]  can  be  generalized, 
viz., 


^abs{ki,Ak  ) 


2567r2q  a3  2 
2kfkf{Ak2)3  ’ 


(4) 


where  a  is  the  fine-structure  constant, 


"*  =  El>  {  E  g/yl) (2£  l,ZLi)+-wl(tLt'L']  WMtt}  ■  V 

In  the  above  l>  is  the  greater  of  l  and  l'  and  W  is  the  Racah  coefficient,  and 

<s  =  So  pv' MJr’  (7) 

where  r_1Pg(r)  is  the  radial  part  of  the  continuum  function  of  the  free  electron  and  V(r)  is 
the  potential  seen  by  the  free  electron.  If  the  dependency  of  s  on  (L,I/,S)  is  suppressed, 

the  quantity  inside  the  curly  bracked  in  Eq.  (6)  reduces  to  M2^,  so  that  Eq.  (4)  reverts  back 
to  Eq.  (6)  of  Ref.  [2], 

The  (mean)  absorption  coefficient  is  obtained  by  averaging  the  cross  sections  over  the 
Maxwellian  distribution  function  f(v,T)  corresponding  to  the  electron  temperature  T,  viz. 


k(A,T)  =  f  Cabsikiy  Ak2)f(vi,  T)dvi, 
Jo 


where 


f(v,T)dv 


4  /  m  \3/2 

V*  [2kTJ  6XP 


(■ 


mv 


•) 


v2dv, 


with 


so  that 


2  kTJ 

=  198.O03/2exp(-31.3296  *:?)&,•  d(fc?), 
0  =  5040/T(iff), 


_  5.05  X  1Q~39Q3/2  r  erp(-31.3296t?)  ;  , 

1  ’  J  (Afc2)3  L  hk, 

The  wavelength  A  and  Ak2  are  related  as 


(8) 

(9) 

(10) 

(11) 

(12) 


A(A)  =  911.8/Afc2(a.u.). 


(13) 


The  cross  sections  in  Eq.(4)  are  expressible  in  units  of  cm5[l]  after  a  substitution  of 
a0  =  0.529xl0_8cm,  so  that  the  absorption  coefficients  in  Eq.  (12)  are  also  in  units  of 
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cm5.  The  absorption  cross  sections  and  coefficients  are  frequently  multiplied  by  the  electron 
density  n  that  would  give  rise  to  a  pressure  of  1  dyne/cm2  at  temperature  T, 


n(cm“3)  =  0.724  x  10l6/T(K) 


(14) 


per  dyne/cm2  [1]  so  that  the  cross  sections  and  coefficients  are  now  expressible  in  cm4/ dyne. 
The  photons  involved  in  the  absorption  also  induce  emission  so  that  the  absorption  coef¬ 
ficients  are  sometimes  multiplied  by  a  factor  [l-exp(-hi//kT)]  [1].  In  this  paper  we  do  not 
include  this  factor  unless  specifically  stated  otherwise. 

The  expression  of  the  emission  cross  sections  has  been  given  in  Refs. [4, 5].  However,  a 
more  pertinent  quantity  in  applications  is  the  emissivity  in  units  of  W-cm3/ ^m-steradian 
[5],  which  we  found  to  be 


J(A,T)  = 


0.794814  x  10 


-31 


03/2  rco  exp(  31.3296  fc,2)^^  (15) 

J  Vo  kj 


[A(/xm)]2 

where  v„  is  the  threshold  velocity  of  the  free  electron  for  emission,  i.e., 


1 

—mv„ 

2  ° 


2  =  hu 


(16) 


and  the  symbol  A(/zm)  indicates  that  A  is  to  be  expressed  in  fim.  The  emissivity  in  this  work 
was  computed  by  Eq.  (15).  Geltman  [5]  also  gave  the  relationship 

1.19  x  104  /  1.439  x  104 


rM  ^  _  1.19  xlO4  (  1.439  x  10^  M  ^ 
“  [A(^m)]«  6XPV  A (jtm)T  ) 

which  can  be  used  to  check  the  consistency  between  Eqs.(12)  and  (15). 


(17) 


2.2  Free  Electron  Function  Py(r) 

The  reduced  radial  function  P*/(r)  is  obtained  as  the  solution  of  the  integro-differential 
equation 

[£  - 

=  E  WifXMVn.i'(r)  +  £feC„.,P„.,(r),  (18) 

n'l'  A  n'V 
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where  V^,  and  are  the  Coulomb  potential  due  to  the  unperturbed  atomic  orbitals  and 
the  polarization  correction  (to  be  discussed  in  Sec.  3)  respectively.  The  first  term  on  the 
RHS  is  the  exchange  interaction  term  and  the  second  term  is  to  ensure  the  orthogonality 
with  the  bound  atomic  orbitals  Vnt.  The  dependency  on  the  angular  momenta  L  and  S  is 
indicated  as  superscripts. 

With  regard  to  the  present  problem  of  0~ [ls22s22p4(3P)kt 2,4 L],  the  unperturbed  orbitals 
of  the  oxygen  atom  are 

#*.(*)  =  Y^typ*  w,  («) 

and  the  Coulomb  potential  is 

=  2-f  -  -  +  2y.(l«.l'W  +  2y.(2»,2j|r) 

+  Ay0(2p,2p\r)  +  ^v£y2(2p,2p|r)}  ,  (20) 

where  Z  is  the  nuclear  charge,  and 

y\(nl,  n'l^r)  =  r'^1  f  Vni{t)Vn>v(t)txdt 

Jo 

+  TX  J^°  Vru{t)'Pn't'{t)t~X~1dt,  (21) 

and  v£  is  equal  to  l/(2i+Z),  -  1,  and  (l+l)/(2l  -  1)  for  L=/+l,  1,1-1  respectively. 

The  exchange  terms  are 

yi  =  2{—  .  yr(l5,H|r)7?i,(r) 

n'l'X  Z4  +  1 

~  ^Y^(25>Wlr)7?2,(r)  +  ^1yi+1(2p,k£\r)V2p(r) 

+  w^f1ye-i(2pyki\r)V2P(r)},  (22) 

where  yt(nl,  kl\r)  can  be  obtained  from  Eq.(21)  by  substituting  P*f  (t)  for  TW(t)- 
The  coefficients  w%s  are  shown  in  Table  1.  Throughout  the  work,  we  have  used  the  wave 
functions  of  dementi  and  Roetti  [17]  for  the  unperturbed  orbitals.  Finally,  Eq.  (18)  is  to 
be  solved  subject  to  the  boundary  condition, 

lim  Ptfir)  =  sin(kr  -  \lv  +  rjj£)t  (23) 
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where  r]fcf  is  the  phase  shift  due  to  interaction  with  the  atom. 

Determination  of  the  polarization  potential  V^(r)  in  Eq.(18)  will  be  taken  up  in  Sec.  3. 
At  this  point  we  should  note  that  if  we  drop  the  V^(r)  term,  then  Eq.(18)  becomes  a  “one- 
channel”  closing-coupling  equation  obtainable  by  omitting  all  coupling  terms  connecting  the 
different  partial  waves  and  the  various  target-atom  states.  The  polarization  potential  may 
be  viewed  as  a  way  to  account  for  the  effects  of  the  omitted  channels  on  the  free-electron 
wave  function.  This  may  be  contrasted  with  the  work  of  John  and  Williams  [7]  in  which  they 
used  the  results  of  a  three-state  close- coupling  calculation  with  no  explicit  introduction  of 
the  polarization  potential  as  the  inclusion  of  the  coupling  channels  in  their  work  is  another 
way  to  allow  for  the  target  polarization. 


3.  METHOD  OF  POLARIZED  ORBITALS 


In  this  section  we  briefly  recall  the  essential  steps  in  Temkin’s  method  of  polarized  orbitals 
[11],  ultimately  leading  to  the  polarization  potential  V^(r)  in  Eq.  (18).  The  details  are 
found  in  the  works  of  Temkin  and  others  [11-13].  Here  the  distortion  of  the  target  atom 
caused  by  the  perturbation  of  the  free  electron  at  f//+ i  is  treated  in  an  adiabatic,  dipole 
approximation  with  the  proviso  that  the  perturbing  electron  is  outside  the  atomic  electrons. 
Thus  the  perturbation  is  taken  as 


N 


'  3  )  r^+1 

e(ritrN+1)  =  j  J’ 


#i=-i 


for  rN+i  >  ri 
for  rN+i  <  r{. 


(24) 


(25) 


To  construct  the  wave  function  of  0(3P)  we  start  with  the  unperturbed  one-electron  functions 
of  the  oxygen  atom  [17]  as  shown  in  Eq.(19).  Following  Temkin  we  add  a  polarization 
correction  to  each  <fi^m  to  form  a  polarized  orbital  which  depends  parametrically  on  r/v+i  as 


<f>nim(r,rN+1)  =  #*(#)  +  X(f>^m{f,rN+1), 


(26) 
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(27) 


*#£»(*»  **+■ l)  =  e(r,rAT+l)  X 

I'm' 


It  is  convenient  to  define 


—  c(im  *  l  m  )Yl'm'(r)rPrUm-*l'm,{.1'), 


(28) 


where 

c{lm-+l'm!)  =  2(-l)m'- 


47t(2£  + 1) 


T  1/2 


(2P  + 1) 


c(l  m'  -mim]  l'm')c{  1  0  l  0;  l'  0).  (29) 


In  the  above  c(jimij2in2'jm)  is  the  Clebsch-Gordan  coeffient;  thus,  we  see  that  | V  — 1  and 
\m!  —  m|=0,l. 

The  radial  function  Pn/m-^wfa)  is  obtained  as  the  solution  of  the  differential  equation 


S3  -  +  E* 


Pnlm^l'm'ir)  =  X  WL"i"\,l'm'(r)Pntm-.l'm'(r) 
n"i"X 

n"l" 

+  rVni(r).  (30) 


Here  the  Coulomb  (V^^)  and  exchange  (W7)  potentials  are  obtainable  by  subjecting  to 
the  Haxtree-Fock  procedure  the  determinantal  wave  functions  in  which  the  parent  orbitals 
^ntm(^)  3X6  replaced  by  the  polarized  version  of  Eq.  (26).  The  function  PnAn-.r,m'(r)  ap¬ 
proaches  zero  as  r  — >  oo,  and  for  small  r 

lirnP„/m^w(r)  =  o/+1.  (31) 

In  this  work  we  have  obtained  the  Pntm—Vm  functions  associated  with  the  unperturbed 
orbitals  nl= 2s  and  2p.  For  later  reference,  they  are  listed  in  Table  2,  and  are  referred  to 
simply  as  P,-  therein.  From  these  polarized  orbitals  we  obtained  the  dipole  polrizability  of 
5.084  a®  compared  with  the  experimental  value  of  5.2±0.4  [18].  Other  theoretical  values 
range  from  4.63  to  5.41  [19].  In  general  the  polarized  orbitals  P„/m_*'m'  depends  on  m  and 
m'  as  well  as  on  l  and  V  due  to  the  difference  in  exchange  term  in  Eq.  (30)  so  that  only  two 
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pairs  are  identical  among  the  fifteen  functions  in  Table  2;  that  is,  P5  and  P6,  and  P10  an-d 
Pi2.  This  dependency  on  m  and  m'  manifests  itself,  for  example,  in  the  calculation  of  dipole 
polrizability  in  which  a  difference  up  to  20%  is  found  in  the  matrix  elements. 

When  the  target  wave  functions  are  described  by  the  unperturbed  orbitals  (j>°,  the 
Coulomb  potential  Vcoul(T)  seen  by  the  free  electron  (coordinate  f),  consists  of  constituent 
terms  due  to  all  occupied  n^-orbitals  such  as 

vMr)  =  (32) 

As  a  consequence  of  the  polarization  of  the  n/^-orbital,  augmentation  of  the  target  orbitals  by 
the  polarization  correction  A <j>^m  in  accordance  with  Eqs.  (26)  and  (27)  produces  additional 
(polarization)  terms,  to  the  first  order  in  A,  of  the  form 

=  /  iS.(r)n.(OP-  '“'r'tfWF'WU. (33) 

The  polarization  potentials  associated  with  the  Coulomb  interaction  which  can  be  derived 
variationally  [11]  are  composed  of  a  series  of  such  integrals,  and  can  be  cast  into  the  general 
form 

t&M  =  -^E  [  ’Pi(t)PmWtdt-  (34) 

The  notation  used  here  requires  some  explanation.  We  have  a  summation  over  j  which  ranges 
from  1  to  15  for  L=£+l,  1  to  21  for  L =£,  and  1  to  31  for  L=£l  as  shown  in  Table  3.  For 
each  value  of  j,  we  also  list  Cy,  u(j),  Vj,  and  B(j)  in  Table  3  for  the  three  possible  choices 
ofL=*+l,  t,  and  t- 1.  To  identify  P s{j)  let  us  take  an  example  of,  say,  j=5  with  L=/+l.  The 
associated  B(j)  is  7  according  to  Table  3.  We  then  go  to  the  entry  of  i=7  in  Table  2  and 
find  that  Pb(j=5)  is  P2P,i-2,2.  From  Table  4,  we  find  Fu^  for  each  j  in  a  similar  way.  Vj 
is  the  reduced  radial  function  of  the  unperturbed  target  orbitals  as  defined  in  Eq.(19)  with 
Vj  =  V2.  for  j<3  and  Vj=V2p  for  j>3. 

The  exchange  potential  is  also  modified  by  the  inclusion  of  polarization.  The  effects  of 
the  polarization  on  the  exchange  term  have  been  examined  by  Henry  [13]  in  his  work  on  the 
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elastic  scattering  of  electron  from  the  oxygen  atom.  The  inclusion  of  this  exchange  term  was 
found  to  change  the  s-wave  (£=0)  cross  sections  by  only  10%.  For  the  higher  partial  waves 
(l  >0),  the  difference  is  even  smaller  (4%)  so  that  the  influence  of  the  polarization  of  the 
orbitals  via  exchange  is  seen  to  be  rather  marginal.  For  this  reason,  we  do  not  consider  in 
this  paper  the  polarization  correction  on  the  electron  exchange. 


4.  NUMERICAL  PROCEDURES 


The  integro-differential  equations  of  Eqs.  (18)  and  (30)  are  solved  non-iteratively  by  using 
the  Lagrange- multiplier  technique  as  we  have  done  previously  [20].  The  inhomogeneity  terms, 
which  arise  due  to  the  electron  exchange,  orthogonality,  and  perturbation,  are  set  to  zero 
after  r=13.5  a0.  This  corresponds  to  the  r  value  at  which  2p  function  has  decreased  to  10-5 
of  its  peak  value.  At  this  point  the  Lagrange  multipliers  are  determined,  and  the  complete 
solution  is  constructed  as  a  linear  combination  of  the  paritcular  integrals  and  homogeneous 
solution.  Beyond  r=13.5  a0  Eqs.  (18)  and  (30)  are  considered  to  be  homogeneous  differential 
equations. 

For  the  free-electron  functions,  the  phase  shift  t)  is  determined  at  various  points  50  to 
250ao  apart  starting  at  450ao,  until  two  successive  values  agree  within  0.01%.  To  solve  for 
the  polarization  correction  functions  of  Eq.  (30),  the  numerical  solution  is  first  initiated 
with  a  certain  arbitrary  value  of  a  in  Eq.  (31).  The  resulting  solution  obtained  by  numerical 
integration  in  general  would  diverge  at  large  r,  say,  tending  to  -foo  as  r— *  oo.  The  value  of  a 
used  here  is  denoted  by  a+.  We  then  change  the  value  of  a  until  the  resulting  solution  tends 
to  -oo  at  large  r,  and  this  new  choice  of  a  is  denoted  by  a_.  By  successive  interpolation 
between  o+  and  a_,  we  develop  a  numerical  solution  for  Pn/m-^'m'(r)  that  monotonically 
decreases  up  to  r=20ao. 
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The  dipole  matrix  elements  are  computed  by  the  velocity  form  instead  of  the  acceleration 
form  shown  in  Eq.  (7).  The  length  and  acceleration  forms  may  be  derived  from  and  is 
equivalent  to  the  velocity  form  if  the  wave  functions  are  exact  eigenfunctions  of  the  atomic 
Hamiltonian.  Using  a  short-hand  notation  of  7  =  (LSk£),  we  have 


I  *,(r) 

Jo 


'  d  (  <>(<'  -  <V 

dr  r 


Py(r)dr, 


(35) 


where  £>  is  the  greater  of  i  and  l' .  The  limits  of  integral  (0,oo)  is  divided  into  two  regions 
(0,R)  and  (R,oo)  where  R  is  the  value  of  r  at  which  the  asymptotic  form  of  Eq.  (22)  sets  in. 
The  integration  in  the  interval  (0,R)  is  carried  out  numerically.  For  the  interval  (R,oo),  the 
integrals  are  expressed  in  closed  forms  analogous  to  those  given  in  Ref.  [21]. 

Using  our  computational  procedures,  we  obtain  virtually  identical  results  for  various  en¬ 
tries  in  Tables  3  and  5  of  Ref.  [2].  In  these  calculations  the  difference  between  the  velocity 
and  acceleration  forms  was  typically  1.5%  or  less. 


5.  RESULTS 

In  Eq.  (6)  we  have  Li  =  Sj=l  for  the  present  problem  (3P  ground  state  of  O  atom), 
and  S  takes  on  1/2  and  3/2.  We  have  also  set  the  limit  of  ( 1,1' )  <4,  and  (L,  I/)  to  all 
values  consistent  with  {1,1').  For  a  given  value  of  the  wavelength  A  (equivalently,  of  Ak2), 
we  have  computed  the  cross  sections  at  200  different  values  of  incident-electron  energy  from 
0  to  1.2  a.u.  as  shown  in  Table  5.  The  process  of  averaging  over  the  Maxwell-Boltzmann 
distribution  in  Eqs.  (12)  and  (15)  was  carried  out  by  interpolating  between  these  200  points. 
The  resulting  continuous  absorption  coefficients  in  Eq.  (12)  are  shown  in  Table  6  and 
emissivity  coefficients  of  Eq.  (15)  in  Table  7. 

Geltman  [5]  made  an  extensive  calculation  of  the  free-free  absorption  coefficients  for  a 
number  of  neutral-atom  systems  including  the  oxygen  atom.  A  Hartree-Fock-Slater-type 
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potential  was  used  to  describe  the  Coulomb  and  exchange  interaction  of  the  free  electron 
with  the  neutral  atom,  and  an  empirical  polarization  potential  was  added.  The  matrix 
elements  are  then  evaluated  in  the  acceleration  form.  At  10,000K,  Geltman’s  values  axe  25- 
30%  smaller  than  our  values  as  shown  in  Table  8.  A  more  detailed  comparison  of  Table  4 
of  Ref.  [5]  and  Table  6  of  the  present  report  shows  that  the  coefficients  in  Ref.  [5]  are  smaller 
than  the  present  values  at  the  lower  electron  temperatures  and  longer  wavelengths  of  the 
absorbed  photon.  At  higher  temperatures  the  difference  becomes  smaller  and  at  20,000K  the 
two  sets  are  nearly  identical.  John  [6]  derived  the  free-free  absorption  coefficients  involving 
various  atoms  and  molecules  from  the  momentum-loss  cross  sections 

by  means  of  the  formula  suggested  by  Dalgarno  and  Lane  [8].  The  data  of  John  for  oxygen 

[6]  are  reported  to  be  in  good  agreement  with  those  of  Geltman  [5],  particularly  for  A  >  1.0/x. 

A  new  and  refined  theoretical  treatment  was  subsequently  offered  by  John  and  Williams 

[7]  in  which  they  applied  the  multichannel  theory  [9]  to  obtain  the  free-free  absorption 
coefficients  for  the  e-0  scattering  system.  For  the  continuum  orbitals  they  used  the  scattering 
data  from  the  three-state  (3P,  1D,  *S  of  the  ground  configuration)  close- coupling  calculation 
of  Saraph  [10].  The  bound  states  were  described  by  the  three- configuration  wave  functions. 
In  Table  8,  we  compare  their  absorption  coefficients  at  9,700K,  computed  by  the  velocity 
form,  with  ours  at  10,000K.  In  making  this  comparison  two  points  should  be  mentioned. 
First,  we  estimate  our  absorption  coefficients  at  9,700K  to  be  only  about  3%  smaller  than  the 
ones  at  10,000K  shown  in  Table  8.  Secondly,  the  coefficients  of  John  and  Williams  contain 
the  contributions  from  cases  other  than  the  oxygen  atom  being  in  the  ground  3P  state  before 
and  after  the  collision.  However,  the  3P-3P  contribution  is  expected  to  occupy  greater  than 
90%  of  the  total  at  9,700K,  since  its  percentage  is  90%  at  12,600K  and  99%  at  5,040K  (See 
Table  3  of  Ref.  [7]).  Thus,  even  after  accounting  for  the  other  contributions,  the  absorption 
coefficients  of  this  work  and  of  Ref.  [7]  differ  by  no  more  than  10%  at  10,000K.  In  Table  5 
of  Ref.  [7],  John  and  Williams  also  give  the  absorption  coefficients  that  include  the  factor 
for  the  stimulated  emission.  In  Table  8  of  this  report  we  include  their  values  at  5,040K, 
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adjusted  so  as  to  exclude  the  stimulated-emission  factor  and  to  be  expressible  in  units  of  cm5. 
These  values  are  seen  to  be  about  10%  smaller  than  our  interpolated  counterparts  at  5,000K. 
Thus,  the  agreement  is  quite  good.  Futhermore,  the  percentage  difference  is  about  the  same 
over  a  wide  range  of  A  and  at  two  different  electron  temperatures.  This  is  encouraging 
and  significant,  because  the  two  formalisms  -  the  close-coupling  method  and  the  method  of 
polarized  orbitals  -  offer  different  approaches  to  allow  for  the  effect  of  target  polarization  in 
a  slow-electron-atom  collision.  This  agreement  may  be  contrasted  with  the  comparison  of 
the  present  calculation  with  that  of  Ref.  [5],  in  which  not  only  greater  differences  are  found, 
but  also  the  differences  depend  on  A  and  T. 

From  a  shock-heated  plasma  experiment  Taylor  and  Caledonia  [14]  obtained  the  absorp¬ 
tion  coefficients,  which  are  shown  in  Table  8  as  “Expt”.  In  a  similar  experiment  Rung 
and  Chang  [15]  found  the  absorption  coefficients  of  7.8  x  10-38  and  3  x  10-36cms  for  A  = 
3.1  and  9.85/4  respectively,  both  at  the  electron  temperature  of  9700K.  Those  experimental 
values  are  about  2  to  3  times  larger  than  the  present  values  and  those  of  Ref.  [7].  Beyond  the 
obvious  difficulties  associated  with  the  experiments  that  are  discussed  in  Refs.  [14]  and  [15], 
we  can  offer  no  clear  explanation  for  the  discrepancy  between  the  experiment  and  theory. 


6.  CONCLUSIONS 

In  an  electron-atom  collision  the  polarization  of  the  target  atom  is  an  important  and 
yet  difficult  feature  to  incorporate  into  the  theoretical  formalism.  The  multi-channel  close¬ 
coupling  method  is  one  of  the  major  theoretical  tools  that  allow  for  the  distortion  of  the 
target  by  the  incident  electron.  In  this  paper  we  approached  the  same  problem  of  the  target 
polarization  with  the  method  of  polarized  orbitals.  We  found  that  this  method  gives  the 
dipole  polarizability  of  5.084  do3  that  is  quite  comparable  to  the  range  of  values  4.63  to  5.41 
by  other  calculations  [19]. 
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With  respect  to  the  free-free  absorption  coefficients  our  values  show  very  good  agreement 
(about  10%)  with  the  results  of  John  and  Williams  [7]  based  on  the  scattering  data  from 
a  three-state  close-coupling  calculation.  This  agreement  is  significant  in  view  of  the  very 
different  nature  between  our  approach  and  that  of  John  and  Williams.  On  the  other  hand 
the  theoretical  absorption  coefficients  are  smaller  by  a  factor  of  about  two  or  three  compared 
to  the  experimental  values,  and  we  are  not  able  to  reconcile  this  discrepancy  at  the  present 
time. 


14 


Table  1.  Values  of  w^s  in  Eq.  (21). 


s 

L 

A  =  *+l 

\  =  l-l 

3 

C'!! 

3(/+l) 

3 1 

2 

'(2/+l)(2/+3) 

(2/— 1)(2/+1) 

1 

*  +  l 

3/(2/+5) 

3/ 

2 

2(2/+l)(2/+3)3 

’  (2/-l)(2/+l) 

X 

l 

3(/-2)  . 

3(/+3) 

2 

2(2/+l)(2/+3) 

2(2/— 1)(2/+1) 

1 

l-l 

3(/+l) 

3(2/— 3)(/+l) 

2 

(2/+l)(2/+3) 

2(2/— 1)J(2/+1) 

Table  2. 

Polarized  corrections  P n^m  that 

are  equivalently 

referred  to  simply  as  P;. 

i 

n  £,m  —* 

l\m! 

i 

n£,m  — » 

1 

2s, 0 

1,1 

9 

2p,l  2,0 

2 

2s, 0 

1,0 

10 

2p,0  2,1 

3 

2s, 0 

1,-1 

11 

2p,0  2,0 

4 

2p,l 

0,0 

12 

2p,0  2,-1 

5 

2p,0 

0,0 

13 

2p,-l  2,0 

6 

2p,-l 

0,0 

14 

2p,-l  2,-1 

7 

2p,l 

2,2 

15 

2p,-l  2,-2 

8 

2P,1 

2,1 

15 


Table  3.  List  of  constituent  parts  that  describe  the 


polarization  potential  V^>/( r)  in  Eq.  (34). 


• 

J 

C; 

u(j) 

Vi 

BO) 

L  =  /  +1 

1 

4/3 

1 

2s 

1 

2 

4/3 

2 

2s 

2 

3 

4/3 

1 

2s 

3 

4 

4/3 

1 

2p 

4 

5 

8/5 

1 

2p 

7 

6 

4/5 

2 

2p 

8 

7 

4/15 

1 

2p 

9 

8 

2/3 

2 

2p 

5 

9 

2/5 

1 

2p 

10 

10 

8/15 

2 

2p 

11 

11 

2/5 

1 

2p 

12 

12 

2/3 

1 

2p 

6 

13 

2/15 

1 

2p 

13 

14 

2/5 

2 

2p 

14 

15 

4/5 

1 

2p 

15 

L  =  / 

1 

4/3 

3 

2s 

1 

2 

4/3 

4 

2s 

2 

16 


Table  3  (2nd  page) 


3 

4/3 

3 

2s 

3 

4 

2/3 

5 

2p 

4 

5 

4/5 

5 

2p 

7 

6 

2/5 

6 

2P 

8 

7 

2/15 

5 

2p 

9 

8 

2/3 

7 

2p 

5 

9 

2/5 

8 

2p 

10 

10 

8/15 

7 

2p 

11 

11 

2/5 

8 

2p 

12 

12 

2/3 

3 

2p 

6 

13 

2/15 

3 

2p 

13 

14 

2/5 

4 

2p 

14 

15 

4/5 

3 

2p 

15 

16 

-2/3 

9 

2p 

5 

17 

-2/5 

9 

2p 

10 

18 

4/15 

9 

2p 

11 

19 

-2/3 

9 

2p 

4 

20 

-2/5 

9 

2p 

8 

21 

4/15 

9 

2p 

9 

L=M 

1 

4/3 

10 

2s 

1 

2 

4/3 

11 

2s 

2 

17 

Table  3  (3rd  page) 


3 

4/3 

10 

2s 

3 

4 

2/3 

12 

2p 

4 

5 

4/5 

12 

2p 

7 

6 

2/5 

13 

2p 

8 

7 

2/15 

12 

2p 

9 

8 

2/3 

14 

2p 

5 

9 

2/5 

15 

2p 

10 

10 

8/15 

14 

2p 

11 

11 

2/5 

15 

2p 

12 

12 

2/3 

16 

2p 

6 

13 

2/15 

16 

2p 

13 

14 

2/5 

17 

2p 

14 

15 

4/5 

16 

2p 

15 

16 

-2/3 

18 

2p 

5 

17 

-2/5 

18 

2p 

10 

18 

4/15 

18 

2p 

11 

19 

-2/3 

18 

2p 

4 

20 

-2/5 

18 

2p 

8 

21 

4/15 

18 

2p 

9 

22 

4/3 

19 

2p 

6 

23 

4/15 

19 

2p 

13 

18 


Table  3  (4th  page) 


24 

4/3 

19 

2p 

4 

25 

4/15 

19 

2p 

9 

26 

2/3 

20 

2p 

6 

27 

-4/15 

20 

2p 

13 

28 

2/5 

20 

2p 

14 

29 

2/3 

20 

2p 

5 

30 

-4/15 

20 

2p 

11 

31 

2/5 

20 

2p 

12 

19 


Table  4.  The  ^-dependent  function  F,(£)  in  the 
polarization  potential  V^(r)  in  Eq.  (34). 


• 

1 

F.M 

i 

F.w 

1 

(M-l) 

ii 

(2P+1-2) 

(2Z+3) 

l(2l-\)(2l+\) 

o 

1 

12 

{4t*+4tz-3P-4l+ti\ 

L 

(2/+3) 

t(2l— 1 )  (21+1)  (2/+3) 

3 

tU+ 2) 

13 

(4/3+Sf2 +9^—15) 

(/+l)(2M-3) 

/(2^-l )  (2/+1)  (2/+3) 

4 

(*+3) 

14 

(4/s+20/*— 13/— 3) 

(/+l)(2/+3) 

/(2/-l)(2/+l)(2/+3) 

5 

/(/+3) 

15 

(4/4+8/s-9/s+3) 

(/+l)(2/+3) 

/(2/-l)(2/+l)(2/+3) 

6 

(1+6) 

16 

(8/<+4/s-8/J+3) 

(l+l)(2i+3) 

/(2/-l)(2/+l)(2/+3) 

7 

1 

17 

(8/s+4/2— 6) 

(*+l) 

/(2/-l)(2/+l)(2/+3) 

8 

/ 

18 

(2*-3) 

(^+1) 

/(2/+l)(2/+3) 

9 

/ 

19 

1 

(/+l)(2/+3) 

(2/+l)(2/+3) 

10 

(2tz-0-l+l) 

20 

(2/-1) 

l(2l-\)(2l+\) 

(2/+l)(2/+3) 

20 


Table  5.  Quadrature  of  the  incident-electron  energy  k2  in  a.u. 
For  each  region  runs  from  kg  in  increments  of  £k2. 


Region 

Id 

£k2 

1 

0.0 

6.25xl0-4 

0.025 

2 

0.025 

1.25xl0-3 

0.050 

3 

0.050 

2.50xl0-3 

0.100 

4 

0.100 

5.0xl0"3 

0.200 

5 

0.200 

l.OxlO-2 

1.200 

21 


Table  6.  Absorption  coefficients  /c(A,T)  of  Eq.  (12)  in  units 
of  cm5.  The  wavelengths  A  are  in  /xm,  and  the  temperature  T 
in  K.  Numbers  inside  the  parentheses  indicate  the  power  of  10. 


X(fim) 

T(K) 

1.0 

9.31(-40) 

1.53(-39) 

2.18(-39) 

2.87(-39) 

4.28(-39) 

6.19(-39) 

2.0 

5.25(-39) 

9.72(-39) 

1.46(-38) 

1.99(-38) 

3.07(-38) 

4.58(-38) 

3.5 

2.30(-38) 

4.63(-38) 

7.20(-38) 

9.98(-38) 

1.56(-37) 

2.37(-37) 

5.0 

6.12(-38) 

1.28(-37) 

2.02(-37) 

2.83(-37) 

4.48(-37) 

6.83(-37) 

7.5 

1.91(-37) 

4.15(-37) 

6.65(-37) 

9.37(-37) 

1.49(-36) 

2.28(-36) 

10.0 

4.35(-37) 

9.66(-37) 

1.56(-36) 

2.20(-36) 

3.50(-36) 

5.36(-36) 

15.0 

1.41(-36) 

3.19(-36) 

5.15(-36) 

7.28(-36) 

1.16(-35) 

1.79(-35) 

20.0 

3.29(-36) 

7.41(-36) 

1.20(-35) 

1.70(-35) 

2.74(-35) 

4.22(-35) 

25.0 

6.30(-36) 

1.43(-35) 

2.32(-35) 

3.30(-35) 

5.32(-35) 

8.21(-35) 

30.0 

1.07(-35) 

2.44(-35) 

3.98(-35) 

5.68(-35) 

9.16(-35) 

1.42(-34) 

22 


Table  7.  Emissivity  coefficients  J(A,T)  of  Eq.  (15)  in  units  of  W  cm3//jm- 
steradian.  The  wavelengths  A  are  in  units  of  fim,  and  the  temperature  T  in  K. 
Numbers  inside  the  parentheses  indicate  the  power  of  10. 


A(/im) 

T(K) 

5000 

10000 

15000 

20000 

30000 

50000 

1.0 

6.22(-37) 

4.32(-36) 

9.94(-36) 

1.66(-35) 

3.10(-35) 

5.24(-35) 

2.0 

4.63(-37) 

1.76(-36) 

3.37(-36) 

5.17(-36) 

8.92(-36) 

1.44(-35) 

3.5 

2.30(-37) 

6.96(-37) 

1.24(-36) 

1.84(-36) 

3.09(-36) 

4.88(-36) 

5.0 

1.31(-37) 

3.66(-37) 

6.37(-37) 

9.35(-37) 

1.55(-36) 

2.43(-36) 

7.5 

6.53(-38) 

1.72(-37) 

2.94(-37) 

4.28(-37)  ! 

7.02(-37) 

1.10(-36) 

10.0 

3.89(-38) 

9.97(-38) 

1.69(-37) 

2.44(-37) 

3.98(-37) 

6.18(-37) 

15.0 

1.88(-38) 

4.57(-38) 

7.60(-38) 

1.09(-37) 

1.77(-37) 

2.75(-37) 

20.0 

1.06(-38) 

2.57(-38) 

4.27(-38) 

6.13(-38) 

9.95(-38) 

1.55(-37) 

25.0 

6.99(-39) 

1.65(-38) 

2.73(-38) 

3.92(-38) 

6.37(-38) 

9.90(-38) 

30.0 

4.77(-39) 

1.14(-38) 

1.89(-38) 

2.72(-38) 

4.42(-38) 

6.88(-38) 

23 


Table  8.  Comparison  of  the  absorption  coefficients  k(A,T)  of  Eq.(14)  in 
cm5  with  other  theoretical  works  (G  and  JW)  and  with  the  experimental 
values  of  Ref.  [14].  Numbers  inside  the  parentheses  indicate  the  power  of  10. 


A(/i) 

2.0 

3.5 

5.0 

10  /x 

Present 
(T  =  10,000K) 

9.72(-39) 

4.63(-38) 

1.28(-37) 

9.66(-37) 

G(Ref.[5]) 

(T  =  10,OOOK) 

7.13(-39) 

8.97(-38) 

6.68(-37) 

JW(Ref.[7]) 

(T  =  9,700K) 

9.6(-39) 

4.67(-38) 

1.31(-37) 

Expt(Ref.[14]) 
(T=9,700K)  ' 

29.0(*39) 

1.0(-38) 

2.6(-37) 

X(fi) 

1.14 

1.52 

2.28 

4.56 

Present 

(T=5,OOOK) 

1.2(-39) 

2.5(-39) 

7.3(-39) 

4.9(-38) 

JW(Ref.[7]) 

(T=5,040K) 

1.09(-39) 

2.28(-39) 

6.60(-39) 

4.38(-38) 
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